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ABSTRACT 

Within the SONYC - Substellar Objects in Nearby Young Clusters - survey, we investigate the 
frequency of free-floating planetary-mass objects (planemos) in the young cluster NGC 1333. Building 
upon our extensive previous work, we present spectra for 12 of the faintest candidates from our deep 
multi-band imaging, plus seven random objects in the same fields, using MOIRCS on Subaru. We 
confirm seven new sources as young very low mass objects (VLMOs), with T e g of 2400-3100K and 
mid-M to early-L spectral types. These objects add to the growing census of VLMOs in NGC1333, 
now totaling 58. Three confirmed objects (one found in this study) have masses below 15Mj up , 
according to evolutionary models, thus are likely planemos. We estimate the total planemo population 
with 5-15 Mj up in NGC1333 is < 8. The mass spectrum in this cluster is well approximated by 
dN/dM oc M~ a , with a single value of a — 0.6 ± 0.1 for M < 0.6 M Q , consistent with other nearby 
star forming regions, and requires a < 0.6 in the planemo domain. Our results in NGC1333, as well 
as findings in several other clusters by ourselves and others, confirm that the star formation process 
extends into the planetary-mass domain, at least down to 6 Mj up . However, given that planemos are 
20-50 times less numerous than stars, their contribution to the object number and mass budget in 
young clusters is negligible. Our findings disagree strongly with the recent claim from a microlensing 
study that free-floating planetary-mass objects are twice as common as stars - if the microlensing result 
is confirmed, those isolated Jupiter-mass objects must have a different origin from brown dwarfs and 
planemos observed in young clusters. 

Subject headings: stars: formation, low-mass, brown dwarfs - planetary systems 
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1. INTRODUCTION 

The frequency of free-floating objects with masses 
below the Deuterium burning limit (< 0.015 M Q ) - 
in the following called 'planemos', short for planetary- 
mass objects - is a subject of debate in the lit- 
erature. Deep surveys in nearby star forming re- 
gions indicate the presen c e of at least some planemos 
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conflicting results regarding their n umbers. On the o ther 
hand, a recent micro-lensing study (jSumi et al.ll20li1) re- 
ports the presence of a large population of Jupiter-mass 
objects in the field ('almost twice as common as stars'), 
which are either without host star or on very wide or- 
bits. Direct detections of field objects in this mass do- 
main have been re ported as well (e.g. iLiu et ail 120111 : 
iLeggett et al.l I2TJ12T) . So far, it is unclear whether the 
mass function declines in the planemo domain or not. 
Young planemos are also useful as benchmark objects 
for testing models for atmospheres and evolution in a 
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new mass and age domain. 

One of the goals of our SONYC survey (short for Sub- 
stellar Objects in Nearby Young Clusters) is to probe the 
planemo domain in several nearby star forming regions, 
using broadband imaging surveys followed by spectro- 
scopic verification of the candidates. One of our targets is 
the young cluster NGC1333, a ~ 1 Myr old compact star 
forming region in the Perseus OB association (see Fig. 
[I]). In this cluster we found a ri ch substellar popu lation, 
with about 30-40 brown dwarfs (|Scholz et al.ll2009L here- 
after SONYC-I). In our most recent paper on this clus- 
ter we also report the discovery of a handful of objects 
with estimated masses below 0.02 M@, one of them at 
roughly 0.006 M Q a nd thus firmly in the planemo regime 
(jScholz et al.H 2012al hereafter SONYC-IV). 

Thus far, the number of planemos in this cluster seems 
low compared with other regions. This, however, is a 
preliminary result, for two reasons: 1) Our spectroscopic 
follow-up was not sufficiently complete in the planemo 
domain. 2) Our optical survey might not be sufficiently 
deep in some of the parts of NGC1333 subject to high 
attenuation due to reddening. Here we aim to resolve 
these issues through additional spectroscopy of faint can- 
didates and a cumulative analysis, thus providing more 
definitive constraints on the planemo population of this 
cluster. 

2. TARGET SELECTION 

The goal of this paper is to probe specifically the 
planemo regime in NGC1333 and to address the poten- 
tial biases outlined in Sect. [TJ To that end, we selected 
three samples of objects for follow-up spectroscopy, in 
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Fig. 1. — Subaru/Suprime-Cam i-band image for our target 
region NGC 1333 in the Perseus s t ar-for ming region. This image is 
identical to Fig. 1 in lScholz et al.l J2012aD . but we have overplotted 
with solid lines the two MOIRCS fields observed for this paper. 
Marked are the cluster radius with a large dashed circle, the two 
well-known objects BD+30 549 (north) and HH 7-11 (south) with 
small ellipses, and the known population of very low mass members 
(Scholz et al. 2009) with small circles. The reflection nebula NGC 
1333 is slightly north of the image center. Coordinates arc J2000. 



total 19 objects, which are listed in Tables Q] and [2] 

1. IZ candidates: Here we started with our primary set 
of 196 candidates identified based on the (i,i-z) colour- 
magnitude diagram (SONYC-I). This diagram is shown 
in Fig. [3J objects for which we previously obtained spec- 
tra are marked with crosses, confirmed very low mass 
objects with squares. A subsample of 22 objects with 
i > 22 has not yet been verified spectroscopically. 20 
of these objects are within 0.25 deg of the cluster cen- 
ter, which is the area where all confirmed very low mass 
members are located. From these 20, we select 6 for 
follow-up spectroscopy, marked with triangles in Fig. [2J 
As can be seen in the figure, these candidates cover i- 
band magnitudes from 22.5 to 25, extending beyond the 
nominal completeness limit of the photometric survey 
(dashed line). Including this new selection, 27 out of 
43 IZ candidates with i > 22 and 9 out of 15 with i > 24 
have been observed spectroscopically. The lowest mass 
confirmed member of the cluster, SONYC-NGC1333-36, 
with an estimated mass of 0.006 M Q appears in this plot 
at i — 23.39 and i — z = 1.96. 

2. JKS candidates: We correlated the deep JK cata- 
logue presented in SONYC-I with the 'HREL' catalogue 
from the Spitzer 'C ores to Disks' (C2D) legacy program 
(|Evans et al.l [20091 . We accepted a source if it has a 
Spitzer counterpart within 1" and uncertainties below 
0.2 mag in the first two IRAC bands at 3.6 and 4.5 /im 
(hereafter II and 12). This yields a catalogue of 824 
objects (called the JKS catalogue). The (J, 11-12) colour- 
magnitude diagram for this sample is plotted in Fig. [3l 
The figure illustrates that this catalogue is significantly 
deeper than the IZ selection: the faintest IZ candidates 
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Fig. 2. — Colour-magnitude diagram for the iz candidates (dots), 
originally identified in SONYC-I. The figure is identical to Fig. 7 
in SONYC-IV, but we have also marked the candidates targeted in 
the new campaign with triangles. Crosses are objects for which we 
have already obtained spectra in previous campaigns. Confirmed 
very low mass objects are marked with squares. The horizontal 
dashed line shows the completeness limit of the photometric survey 
as estimated in SONYC-I. 
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Fig. 3. — Colour-magnitude diagram for the JKS candidates 
(dots), see text for details. Confirmed very low mass objects are 
marked with squares, crosses are objects for which we have already 
obtained spectra in previous campaigns. Objects targeted in this 
new campaign are marked with triangles. The completeness limit 
of the J-band photometry is marked with a dashed line. 

have J ~ 19, whereas the JKS catalogue is complete 
down to J = 20.8. From this catalogue, we only consider 
the subsample with 71 — 72 > 0.3. This cut-off was intro- 
duced to avoid the bulk of the datapoints in the diagram, 
for which a spectroscopic follow-up is not practical; in 
Sect. 15.31 we will discuss the consequences of this cut on 
the completeness of the survey. There are only 17 objects 
without spectroscopic verification with 16 < J < 19, 
II — 12 > 0.3 and (as above) a distance from the clus- 
ter center below 0.25 deg. Out of these 17, we selected 6 
for follow-up spectroscopy, which are marked with trian- 
gles. The aforementioned SONYC-NGC1333-36 appears 
in this plot at J = 18.53 and II - 12 = 0.35. Note that 
spectroscopic verification for objects beyond J = 19 is 
not feasible with the currently available facilities. 

3. Random candidates: After selecting the IZ and JKS 
candidates, we decided to observe two MOIRCS fields 
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(see Sect. [3] and Fig. []}, which cover a large fraction 
of the samples. To prepare the spectroscopy run, we 
took a new set of deep K-band images for these two 
fields. From these images, we selected 7 more sources 
for spectroscopy, for which we had no prior information 
about their colours. The only constraint for these ran- 
dom sources was the outline of the slit mask. 

We note that most of the selected objects are 
listed as candidates from pho tometric su rv eys i n 
the literat u re, i n papers by Aspin et al.l (11994ft; 
lLada et all (11996ft; Iwilking et al.l (12004ft; lOasa et al.l 



( 2008lT lGutermuth et al.l (12008ft . 

Two caveats need to be pointed out with regard to the 
selection of candidates from colour-magnitude diagrams, 
as done here for the IZ and JKS candidates. 

First, the completeness of the survey depends strongly 
on the definition of the colour and magnitude criteria 
used to select the candidates for follow-up spectroscopy. 
For the IZ sample, we defined the selection box based 
on the position of previously identified brown dwarfs and 
extend it at the faint end - we are simply probing the ex- 
tension of the colour space occupied by the known brown 
dwarfs towards fainter magnitudes. A similar strategy 
was used for the JKS sample, but as argued above, we 
need to introduce an additional colour cut (see also the 
discussion Sect. 15.31) . In the two samples, the colour cuts 
are independent of evolutionary models. 

Second, the photometric uncertainties may cause us to 
miss some brown dwarfs. For the iz selection, this is not 
a major problem: The median error (excluding system- 
atics) for i > 22 is 0.05 mag, and 95% of the errors are 
below 0.1 mag. The width of the selection box (lmag 
in (i-z) colour) is significantly larger the errors. In ad- 
dition, most of the confirmed objects are 0.1-0.2 mag to 
the right side of the left boundary of the box. Assuming 
that these objects show a reliable representation of the 
cluster isochrone in NGC1333, even large uncertainties 
are unlikely to move the datapoints of substellar mem- 
bers out of our selection box. Note that the right side of 
our selection criterion was only formal, no objects in our 
survey had colours redder than this boundary. 

In the JKS sample, the median error in the (11-12) 
colour is 0.1 mag for 16 < J < 19. For 18 < J < 19, 
i.e. at the faint end of our candidate selection, the errors 
are < 0.22 mag. These errors do not depend significantly 
on the (11-12) colour. The expected (11-12) colours of 
brown dwarfs and planemos are close to the bulk of the 
datapoints in this diagram (see Sect. 15. 3[) , i.e. we may 
lose objects on the left side of our colour cut due to large 
photometric errors. 

3. OBSERVATIONS AND DATA REDUCTION 

New spectra for candidate planemo members of 
NGC1333 were taken using the near- infrared spectro- 
graph MOIRCS at the Subaru telescope in multi-object 
(MOS) mode. Since our targets are faint and require 
long integration times, MOS is the most practical way of 
conducting a spectroscopic survey. After creating the IZ 
and JKS candidate samples, we looked for the optimum 
way to obtain spectra for as many as possible of these 
objects. We found that two MOIRCS fields arranged as 
shown in Fig. [T] can cover 12 of these candidates (6 IZ 
and 6 JKS, as discussed in Sect. [2]). 

K-band pre-imaging was carried out for these fields in 



September 2011, with excellent seeing (0.3") and under 
good conditions. The pre-images were reduced and cali- 
brated using 2MASS photometry. Aperture magnitudes 
were measured for all our candidates, to complement and 
verify the available photometry. These pre-images are 
used to define the masks for the spectroscopy run, which 
was carried out in the night of Nov 14-15. For most of 
this night, we observed under stable conditions and with 
seeing in the range of 0.5" . At 2am local time, about 2 
hours before the end of the visibility for our target region, 
the humidity increased above 80% forcing us to close the 
dome. In total, we spent 3:20 h telescope time on mask 
1 and 3:10 h on mask 2. At the beginning of the night 
and between mask 1 and mask 2 we observed 3 standard 
stars with spectral type around AO intended to be used 
for telluric calibration. 

We used our two MOS masks, the grism HK500 with a 
nominal resolution of 600-800 and a wavelength coverage 
from 1.3 to 2.4/zra, an order sorting filter, and slits with 
a length of 11.7" and a width of 1.17" . For each mask we 
obtained a series of individual exposures of at most 5 min. 
Due to an incident with the spillage of cooling liquid ear- 
lier in the year (July 2 2011), the autoguider system was 
not operational at the time of the observations. There- 
fore, we checked the position of the mask alignment stars 
between individual exposures and re-aligned the mask if 
the offset was more than the typical seeing in the direc- 
tion perpendicular to the slit. If the offsets were large, 
we reduced the integration time to be able to track the 
sources reliably. The mask position was typically stable 
over several minutes, the integration times for individ- 
ual exposures were 240 sec for mask 1 and 80-300 sec for 
mask 2. The total on-source time was 88 min for mask 1 
and 60 min for mask 2. Between consecutive exposures, 
the telescope was nodded along the slits by 2" . 

We reduced the MOIRCS spectra using the recipes de- 
scribed in detail in SONYC-I. This includes the a) sub- 
traction of consecutive exposures that are offset along 
the slit, b) division by a normalised domcflat calculated 
from lamp-on and lamp-off exposures, and c) the stack- 
ing of the individual frames for each mask. For the last 
step, we measured the position of one of the sources in 
the slit for each individual frame and shifted the frames 
accordingly, to account for small offsets. 

Spectra were extracted using apall in IRAF. This in- 
cludes aperture definition, background subtraction, and 
trace fitting. We opted for a full aperture integration 
across the profile (as opposed to optimum extraction), 
because many of our sources are faint and lack a well- 
defined profile. The same procedure was carried out for 
the standard stars. The wavelength solution was deter- 
mined from ~ 10 telluric lines in the raw spectra. For 
each target, the two resulting spectra, corresponding to 
the two positions in the slit, were coadded. The final 
spectra are re-binned to 40 A per wavelength element. 

The three standards HIP6855, HIP9196, and HIP21115 
have been observed at airmasses of 1.46, 1.64 and 1.10, 
while the average airmass for the science exposures was 
1.30 for mask 1 and 1.06 for mask 2. The spectrum for 
HIP6855 has a significantly lower flux level in the K- 
band (by 10-20%) than the two other stars and was not 
used. The scaled spectra for HIP9196 and HIP21115, 
however, are very similar with small differences in the 
range of 1-2%, well-explained by the differences in air- 
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mass. This is consistent with an extinction coefficient 
difference (mag/airmass) in the range of 0.01 between 
H- and K-band, which is plausible for Mauna Kea. 

Because HIP21115 was observed between the two 
masks and thus should give the best approximation to 
the conditions during the science exposures, we used its 
spectrum for the calibration. All science spectra were di- 
vided by the quotient between the telluric standa rd spec- 
trum and the tabulated spectrum of an AO star ()Picklesl 
[19981) . 7 

4. SPECTRAL ANALYSIS 

The 19 available spectra were analysed using the 
recipes described in detail in SONYC-IV. In the follow- 
ing, we give a brief summary of the four steps and de- 
scribe the results. 

1. Selection of young very low mass objects: 7 
spectra show the characteristic spectral signature of 
young very low mass objects. Specifically, they exhibit 
the sharp 'H-band peak' at 1.68 /im caused by broad 
water absorption fea tures at 1.3-1.5 and 1.75-2.05 fim 
(iCushing et al.l 120051). which is an indication for low- 
gravitv (|Brandeker et al.ll2~006D . These 7 sources are con- 
sidered to be very low mass objects. Given the compact 
nature o f NGC1333 and the low space densities for these 
objects (ICaballero et al.ll2008l ). contamination by fore- 
ground or background objects is negligible. Thus, these 
7 objects are most likely young members of NGC1333. 
Their spectra are shown in Fig. [4j their properties listed 
in Table [T] The remaining 12 spectra are either feature- 
less or show only little structure in the H-band. They 
are listed in Table [2J 

2. Spectral types: For the 7 objects from Table [1] we 
measured spectral types using the H-peak index (here- 
after HPI) introduced in SONYC-IV, after dereddening 
the spectra using the spectroscopic Ay (see below). As 
discussed in the previous paper, other near-infrared in- 
dices suggested in the literature are not suitable for our 
type of data. The HPI calculates the flux ratio between 
the wavelength intervals 1.68 and 1.50 [im and is usable 
for >M6 spectral types. Using the calibration given in 
SONYC-IV, all 7 objects identified here have spectral 
types of M7 or later (see Table [T]), the latest type is ~L1 
(albeit with a noisy spectrum). The typical uncertainty 
is one subtype. 

3. Effective temperatures: We calculate a grid of model 
spectra for effective temperatures from 1500 to 3900 K 
(in steps of 100 K) and extinctions Ay from to 30 mag. 
This g rid is based on the DUSTY models bv lAllard et al.l 
(2001). For all objects with any indication of struc- 
ture in the H-band, we searched for the best match- 
ing model spectrum. This was done in the following 
way: For each model we subtracted observed spectrum 
from model spectrum, squared this quantity, divided by 
the model spectrum, and averaged the residual over the 
wavelength domain 1.38-1.78 and 2.10-2.32 /im. We se- 
lected the model spectrum for which this quantity is min- 
imal, and adopted its parameters as best fit (see Table 
[T]). The typical uncertainty is ±200 K in T e g and ±1 mag 

7 Although HIP21115 is a B9V star, we used the tabulated AO 
spectrum because the Pickles et al. spectrum for spectral type 
AO has much better resolution. The difference in the flux levels 
between AO and B9 do not matter in the H- and K-band for our 
purposes. 
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Fig. 4. — Near-infrared spectra of the 7 newly identified very low 
mass members of NGC1333 (solid lines). The dashed lines show 
the best-fitting, reddened model spectra with the parameters as 
listed in Table [T] The spectra are offset on the y-axis by constants 
for clarity. 



in Ay. Whenever possible we also calculated Ay from 
the J- and K-band photometry. The results are given 
in Table [TJ The 7 very low mass objects have effective 
temperatures between 2400 and 3100 K. The remaining 
sources with recognisable structure in the H-band have 
temperature around or above 3500 K (Table \2§ ■ The ob- 
ject SONYC-NGC1333-42, the coolest source in our sam- 
ple, has a noisy spectrum; therefore we manually chose 
the best-fitting model. The best-fitting reddened models 
are overplotted in Fig. 01 

Two of the objects in Table Q] were selected randomly. 
One of them (SONYC-NGC1333-45) turns out to be on 
our initial list of IZ candidates but much brighter than 
the cutoff used for this paper [i = 18.609 vs. i = 22) . 
It was previously confirmed by IWinston et al.l (|2009j ). 
They find a spectral type of M8 and a temperature of 
2700 K; our analysis yields consistent results (M8 and 
2900 K). The second random source, SONYC-NGC1333- 
44, has been listed before as a pho tometric candidate 
(|Aspin et all 11994 lOasa et ail 120081) . It was not found 
in our IZ survey because it sits very close to a blooming 
spike caused by a bright star. 

4- Sanity checks: As in our previous papers, we 
additionally compared the spectra to literature spec- 
tra of other brown dwarfs. This is particularly impor- 
tant to make sure that our results are not severely af- 
fected by the often low signal-to-noise ratio in our data. 
For comparison we use the spectra for the three young 
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TABLE 1 

New very low mass members of NGC1333 



ID 


a(J2000) 


<5(J2000) 


J (mag) 


K (mag) 


SpT 


A v (phot) 


A v (spec) 


T e g 


Comments 


39 


03:28:54.96 


+31:18:15.3 


16.881 


13.155 


M8.5 


15 


17 


2900 


IZ 


40 


03:29:02.80 


+31:22:17.3 


16.959 


13.423 


M9.3 


14 


15 


2700 


JKS 


41 


03:29:04.63 


+31:20:28.9 


17.647 


13.892 


M7.6 


15 


18 


3100 


JKS 


42 


03:29:13.85 


+31:18:05.7 


19.252 


17.540 


LI 


4 


1 


2400 


IZ a 


43 


03:29:21.94 


+31:18:29.2 


18.728 


17.128 


M8.1 


3 


5 


2900 


IZ 


44 


03:28:57.04 


+31:16:48.7 




15.59 


M8.6 




4 


2800 


random 


45 


03:29:13.04 


+31:17:38.3 


15.23 


14.16 


M8.1 





2 


2900 


random 



Fit manually adjusted to account for poor spectral quality 



TABLE 2 

Objects ruled out as very low mass 
members of ngc1333 



a(J2000) <5(J2000) 



Comments 



03:29: 
03:28: 
03:29: 
03:29: 
03:28: 
03:29: 
03 29 
03:29: 
03:28: 
03:29: 
03:29: 
03:29: 



04.31 
59.33 
01.88 
02.69 
58.42 
05.67 
13.05 
00.02 
56.31 
06.45 
16.26 
21.88 



+31:19:06.4 
+31 16 31.5 
+31:16:53.4 
+31 19 05.6 
+31:22:17.6 
+31:21:33.9 
+31 16 39.8 
+31:17:13.4 
+31:20:35.2 
+31:16:52.0 
+31:18:06.6 
+31:17:40.7 



JKS, T - 3500 
IZ, T > 3500 
JKS, T ~ 3500 
IZ, T ~ 3500 
JKS, featureless 
JKS, T - 3500 
IZ, T > 3500 
random, T > 3500 
random, featureless 
random, featureless 
random, featureless 
random, featureless 



ob jects 2M1207A DHT auB, and 2M1207B, published 
bv iPatience et all ()2012j ) in a detailed comparison with 
model spectra. For all three they find a satisfying match 
with different types of model spectra and consistent tem- 
peratures for a variety of spectral ranges and model spec- 
tra. According to their analysis, these objects have ef- 
fective temperatures of 3100, 2600, and 1600 K, when 
HK spectra are compared with DUSTY models, which 
is comparable to the way we have determined tempera- 
tures. With our method, we derive T e g of 3200, 2800, and 
1700 K for the same data, consistent with the published 
values. 

In Fig. [5] we plot their spectra in comparison with 
three of our newly identified objects in NGC1333 with 
similar temperatures: SONYC-NGC1333-41, -40, and 
-42 with T cff of 3100, 2700 and 2400 K. The SONYC 
spectra were dereddened using the spectroscopic Ay 
given in Tabic Q] From Fig. [5] it seems plausible 
that the two pairs SONYC-NGC1333-41/2M1207A and 
SONYC-NGC1333-40/DHTauB have similar tempera- 
tures. Based on the slope of the H-band peak, it is also 
justified to say that SONYC-NGC1333-42 is probably 
somewhat hotter than 2M1207B. There are some offsets 
in the flux level in the K-band, which may be explained 
by inaccuracies in the extinction and/or excess flux from 
a disk. Overall, the results of our analysis are confirmed. 8 

For the same three new SONYC objects, we compare 
the photometric fluxes from 0.8 to 8 //m, as far as avail- 



8 The study by Patience ct al. (2012) includes three more young 
sources - CTCha, GQLupB, TWA5B. For these objects their 
comparison with models does not yield consistent temperatures 
and the fit with the model spectra is poor. Therefore we do not 
use them for the comparison, but these objects clearly pose a chal- 
lenge for our understanding of substellar spectra. 
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Fig. 5. — Three spectra for newly identified very low mass mem- 
bers of NGC1333 from Table [T] (solid lines) in comparison with 
published spectra for brown dwarfs from Patience ct al. (2012) with 
similar effective temperature (dashed lines). The spectra are offset 
on the y-axis by constants for clarity. 

able, with the model spectra, adopting T c ^ and the pho- 
tometric Ay as listed in Table [TJ see Fig. [5] Again, 
this comparison confirms the results from our analysis. 
For SONYC-NGC1333-40 mid-infrared excess is clearly 
visible and can be attributed to the presence of a disk. 
As pointed out in SONYC-IV, spectroscopy covering the 
wide spectral range from the J- to the M-band promises 
to be the ideal tool to improve the characterisation of 
substellar objects. 

5. DISCUSSION 

With the new spectra presented in this paper we have 
enlarged the number of known very low mass members 
in NGC1333. In SONYC-IV we list 51 confirmed objects 
with spectral type M5 or later and/or effective temper- 
ature of 3200 K or lower (41 previously known and 10 
new); this new campaign adds 7 new objects to this cen- 
sus. Based on the currently available isochrones for this 
mass domain, most of these 7 are cool enough to be con- 
sidered to be brown dwarfs. In the following we will 
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Fig. 6. — Photometric spectral flux distribution for three newly identified very low mass members of NGC1333 in comparison with 
DUSTY model spectra for T e g as given in Table [l] The model spectra have been dereddened using the photometric Ay listed in Table [T] 
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Fig. 7. — Distribution of spectral types for the very low mass 
population in NGC1333, including the previously confirmed ob- 
jects with >M5 listed in SONYC-IV, the new objects found in 
SONYC-IV and this p aper, and the earlier type objects classified in 
I Winston et al. (2009). Note that Winston et al. additionally clas- 
sify 12 objects with spectral types earlier than M0. Approximate 
mass limits according to DUSTY/COND isochrones are given. 

discuss the implications of the new results for the sub- 
stellar population in NGC1333, with particular focus on 
the plancmo domain. 

5.1. The distribution of spectral types and magnitudes 

We begin by showing the spectral type distribution for 
all spectroscopically classified M and L-type objects in 
NGC1333 in Fig. [7] This includes the very low mass 
objects listed in SONYC-IV, new objects identified in 
the same paper a nd this new stu d y, plu s earlier type ob- 
jects classified in IWinston et al] (|2009t ). Their spectro- 
scopic survey is essentially complete down to the brown 
dwarf/star boundary (for Ay < 10), and thus ideally 
complements our own studies, which are not complete 
above the substellar boundary. 

For the stellar regime, the distribution of spectral types 
in NGC1333 resembles th ose found for other regions 
(Chal, IC348, Taurus, see lLuhmanll2007l ). with a peak 
around M4 and a increase in the frequency between M0 
and M5 by about a factor of 3. For later spectral types, 
we observe a second peak around M7, which is not seen 
in these three other regions. 

Another relevant feature in the context of this paper 
is the s harp drop of the nu mbers at M9. Based on the 
CON E) (jBaraffe et al.ll2003h and DUSTY (jChabrier et al.1 
I2000T ) isochrones, M9 corresponds to a mass around the 
Deuterium burning limit. Only two objects (SONYC- 
NGC1333-36 at ~L3 and SONYC-NGC1333-42 at LI) 
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Fig. 8. — Distribution of absolute J-band magnitu des for 
the Class II in NGC1333 from IGutermuth et al.1 12003) for all 
sources with Ay < 10. Approximate mass limits according to 
DUSTY/COND isochrones are given. 

have later spectral types. Thus, the number of observed 
planemos in this cluster is very low. This feature will be 
further discussed in the following subsections. 

As an independent check of the distribution shown in 
Fig. we examine the histogram of absolute J-band 
magnitudes for the pop ulation of Class II sour ces in 
NGC1333, as identified bv lGutermuth et all (|2008l ) based 
on Spitzer data. Of their 137 objects, 94 have photom- 
etry in 2MASS. We calculate Ay from J- and K-band 
photometry, using the same prescription as in SONYC- 
IV, deredden the J-band photometry and subtract the 
dis tance modulus of 7. 4 (assuming a distance of 300 pc, 
see lBelikov et al.1 1200 2). We exclude the most reddened 
objects and consider only the 74 with Ay < 10. In this 
extinction regime the distribution of J-band magnitudes 
does not show a dependence on Ay. 

Fig. [8] shows again a peak slightly above the substellar 
limit, which corresponds to the peak around M4-M5 in 
the spectral type distribution. A second peak is visible 
in the substellar regime (for Mj of 7-8), but it is much 
less pronounced than the peak at M7 seen in Fig. [7] In 
fact, within the statistical errors the magnitude distribu- 
tion is consistent with a broad peak around Mj = 5 and 
a decline at fainter magnitudes, which agrees with pre- 
vious ly published results for Cha-I and IC348 ()Luhmaij 
2007). Thus, the second peak in the spectral type distri- 
bution could be spurious. One possible explanation is a 
gap i n completeness betwe en the spectroscopic samples 
from IWinston et all f|2009T ) and ours, causing us to miss 
objects at ^M6. 
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Fig. 9. — HRD for the very low mass population in NGC1333 
(crosses - previously known objects in the census from SONYC-IV, 
plusses - new objects found in SONYC-IV, squares - new objects 
identified in this paper) . The circle shows the improved parameters 
for SONYC-NGC1333-36 after fitting the photometric SED, see 
SONYC-IV. Model isochrones (solid - BCAH, dotted - DUSTY, 
dashed - COND) are shown for 1 (upper) and 5Myr (lower). The 
dash-dotted lines show the limits of our spectroscopic survey for 
Ay = and 5 mag. 

5.2. The Hertzsprung-Russell diagram 

In Fig. |9]we show the HRD for NGC1333, focusing on 
the very low mass objects with spectral types of M5 or 
later. Out of the 41 previously known sources listed in 
SONYC-IV (crosses), the 10 newly identified in SONYC- 
IV (plusses), and the 7 newly identified in this paper 
(squares), 51 have estimated effective temperatures and 
are plotted in the figure. 

Instead of luminosity, we prefer to plot the absolute 
J-band magnitude to avoid additional errors from the 
bolometric correction. The absolute J-band magnitudes 
are calculated in the same way as in Sect. 15.11 includ- 
ing dereddening of the J-band data 9 and subtraction of 
the distance modulus. The typical errorbar in the ab- 
solute J-band magnitude is ±0.5 mag, which combines 
the errors in Ay (±lmag), distance (±20pc), and pho- 
tometry (up to ±0.1 mag). The uncertainty in effective 
temperature, as outlined above, is ±200 K. The approx- 
imate mass limits according to the DUSTY and COND 
isochrones, assuming an age of 1 Myr, are given as well. 

The general trend in this diagram is in good agreement 
with 3 sets of isochro nes from the Lyon g roup: DUSTY, 
COND, and BCAH (|Baraffe et all [1998) . Typical for 
HRDs of star forming regions, th e objects show a larg e 
spread due to various effects (e.g. iLittlefair et alj|201lj ). 
However, only 7 objects (14%) are inconsistent with the 
1-5 Myr isochrones within the errorbars. Some of them 
could be affected by ongoing accr etion or strong ma, 
netic activity (overluminosity) - see lScholz et al.l ((2012 
for a discussion - or grey extinction by an edge-on disk 
(under lumino sity) . Some might experience variability 
(Scholz 2012). In addition, there remains a small chance 
that an outlier is a fore- or background so urce: Perseus 
is known to have a complex substructure (|Enoch et al.l 

9 As in SONYC-IV, we use an intrinsic J-K colour of 1.0, which 
is roughly applicable down to the planemo regime. For comparison, 
the l atest type objects in a Ori (7 objects) have an average J-K of 
1.25 ICaballero et all 12007). A shift in J - K from 1.0 to 1.25 
causes a shift of 1.36 mag in Ay. 



120061: iRidge et al.ll2006D . i.e. there may be young objects 
at different distances. Moreover, in some cases the signal- 
to-noise in our spectra is not sufficient to confirm low 
gravity, i.e. we cannot completely exclude older objects 
in the field. Given these caveats, the overall agreement 
with the isochrones is good. 

Only a very small number of 3 objects have tempera- 
ture below 2500 K - SONYC-NGC1333-31, 36, and 42 - 
(two of them also with spectral type later than M9, see 
Sect. 15. ip and would have masses around or below the 
Deuterium burning limit. The plot demonstrates that 
the frequency of objects per bin in effective temperature 
drops significantly below 2500 K. There are 19 objects 
with 2800 < T cff '< 3100, 13 with 2500 < T ff < 2800, 
but only 3 with 2200 < T cS < 2500. This confirms the 
finding from the spectral type distribution. 

The drop in the frequency of observed objects below 
2500 K and spectral type M9, already found in Sect. 15. 11 
can either be explained by an actual lack of planemos in 
this cluster, or by biases in our photometric survey and 
the selection of objects for follow-up spectroscopy, which 
will be discussed in the following subsection. In addition, 
there are three more caveats: 

1) Extinction might limit the completeness of our sur- 
vey below 2500 K. To test this, we have overplotted with 
dash-dotted lines the J-band limit of our survey (J = 19, 
see Sect. [2]) for Ay — and Ay — 5. We should have 
been able to find objects with T e g > 2200 K and Ay < 5. 
The clear majority of the confirmed very low mass ob- 
jects have Ay < 5. Also, Ay ~ 5 is the median ex- 
tinction for the sample of Class II objects published by 
iGutermuth et alj (|2008f) . as far as their 2MASS photom- 
etry is available. Therefore we are unlikely to miss a sig- 
nificant number of planemos down to 2200 K due to their 
high extinction. For extinctions well below Ay of 5 mag 
we should be able to find sources even cooler than 2200 K 
(as proven by the detection of SONYC-NGC1333-36). 

2) The isochrones might not accurately reflect the lu- 
minosity vs. effective temperature relation in young star 
forming regions. To explain a lack of planemos in Fig. 
O the isochrones would have to turn downwards below 
2500 K, i.e cooler objects need to be significantly fainter 
than expected based on these models. While this cannot 
be definitely excluded at the moment, the data for the 
eclipsing binary brown dwarf in Orion - currently the 
only available empirical te st for these isochron es - does 
not indicate such a trend (|Stassun et al.|[200~7T ) . The lu- 
minosities for its components arc in agreement with the 
predictions from the isochrones, in fact, for the lower 
mass component (M = 0.036 Mq) the isochrone under- 
estimates the luminosity, the reverse of what would be 
needed to explain the lack of planemos in NGC1333. Fur- 
thermore, it is difficult to conceive a significant departure 
from the smooth (almost linear) trend in the isochrones 
in Fig. 31 since this would require an additional and 
strong source of opacity in the J-band. To our knowl- 
edge, this is not born out by the available data, for exam- 
ple, the same isochrones fit the optical and near-infrared 
colours for the substellar p opulation in a Ori qui te well 
down to the lowest masses (|Caballero et alJl200l . 

3) The effective temperatures - mostly determined 
with the routines used in this paper, i.e. by fitting 
DUSTY models to the near-infrared spectra - could be 
systematically unreliable in a way that T c g is overesti- 
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mated for very cool sources. The study bv lPatience et al.l 
(|2012T ) . already mentioned in Sect. Ill does not in- 
dicate any such trend. Specifically, their comparison 
the DUSTY models yield effective temperatures below 
2000 K when compared with the spectra for 2M1207B, 
CT Cha B and GQ Lup B, three very young objects with 
L spectral types. This seems to work well even when only 
the H- and K-band spectra are considered, which is the 
case for our objects in NGC1333. 

We conclude that these three caveats cannot be used to 
explain the lack of observed objects cooler than 2500 K 
and later than M9. 

5.3. Survey completeness 

In this subsection we will explore the possibility that 
our survey has missed a significant number of planemos 
due to a lack of depth or biases in the colour cuts. 

Our initial candidate selection was based on a (i, i-z) 
colour magnitude diagram (the IZ catalogue, see Sect. 
[2]). From the set of 196 candidates, more than half are 
now covered by follow-up spectroscopy. In particular, 
27 out of 43 IZ candidates with i > 22 and 9 out of 
15 with i > 24 have been verified spectroscopically. As 
illustrated in Fig. [21 the spectroscopy covers the entire 
colour- magnitude space of this selection down to the faint 
end. This candidate sample contains 2 likely planemos 
(SONYC-NGC1333-31 is not found in this survey). Sta- 
tistically, the remaining candidates which have not been 
verified yet are expected to contain one more. 

This selection, however, might have missed planemos, 
because the depth of the IZ survey is not sufficient. As 
argued in SONYC-IV, this is a distinct possibility; while 
the nominal completeness limit is i = 24.8, it might be 
lower in some parts of the cluster with strong background 
emission. 

To overcome the depth limitation of the IZ survey, we 
have used the JKS catalogue, which contains all objects 
with J- and K-band photometry from our own survey 
and Spitzer photometry in the first two IRAC bands at 
3.6 and 4.5 [im (see Sect. [2j). As illustrated in Fig. [3l 
this survey is much deeper and extends significantly be- 
yond the limiting magnitude for follow-up spectroscopy 
(J = 19). It should be complete down to J = 20.8. 
The near-infrared images do not show strong background 
structure due to the cloud, i.e. the completeness should 
be uniform across the cluster. The JKS catalogue con- 
tains about two thirds of the confirmed very low mass 
objects in NGC1333; most of the ones that are not in 
this catalogue are relatively bright and therefore affected 
by saturation in our deep near-infrared images. This cat- 
alogue contains 2 likely planemos (SONYC-NGC1333-42 
does not have a Spitzer counterpart). 

The colour-magnitude diagram from this catalogue is 
shown in Fig. [31 Our spectroscopic follow-up covers 
a large fraction of this diagram to the J = 19 limit. 
According to the COND and DUSTY 1 Myr isochrones, 
this should correspond to a mass of 0.003-0.004 M for 
A v = and 0.006-0.008 M e for A v = 5. For 15 < 
J < 17 the entire colour range in II — 12 has been 
probed by spectroscopy; it turns out that most of the 
confirmed objects in this magnitude range (11/13) have 
II - 12 > 0.2. As the IRAC colour can be affected by 
excess emission due to a disk, this simply might reflect 
the high disk fraction in NGC1333 (see SONYC-IV). For 
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Fig. 10. — Colour-colour plot for the JKS catalogue (see Sect. 
[2{. Marked are objects for which we have taken spectra (crosses) 
and confirmed very low mass members of NGC1333 (squares). The 
dash-dotted lines mark the colour space where brown dwarfs and 
planemos are expected to be found. The dashed line shows the 
reddening path for the average planemo. 

17 < J < 19 we verified about three quarters of the ob- 
jects with II — 12 > 0.3. Thus, in these colour domains 
the spectroscopic survey is sufficiently complete to rule 
out the existence of more than 1-2 additional planemos. 

To further analyse the colour space from the JKS cat- 
alogue, we show a (J-K, 11-12) colour-colour diagram 
in Fig. 1101 In comparison with other colour combina- 
tions in this catalogue, this plot shows a relatively clear 
distinction between objects confirmed as very low mass 
cluster members and those rejected. All confirmed ob- 
jects are marked with squares, all objects verified with 
spectroscopy with crosses. For this plot, we estimate 
the expected position of planemos based on the pho- 
tometry ^_of_k2iow[n_youn^ In the a Ori 
cluster (jScholz fc Javawardhanal2008| ) the average colour 
of 16 M9-L5-type objects is II - 12 = 0.2, albeit with 
a spread of ±0.2 mag due to photometric uncertainties. 
Among 10 young field dwarfs with L0-L4 spectral type, 
the avera ge colour is II — 12 = 0.17 with a spread of 
0.05 mag (jLuhman et al.ll2~009D . Based on these two sam- 
ples, we adopt II — 12 = 0.2 as typical intrinsic colour of 
planemos. 

The dashed line in Fig. [TU] shows the reddening path 
for Ay = — 15 starting with this intrinsic colour. Shift- 
ing this line by 0.1 mag to the left is in good agreement 
with the left boundary for most of the confirmed very 
low mass objects in NGC1333 (dash-dotted line). This 
line combined with a limit of J — K > 1.0 is adopted as 
the colour space where brown dwarfs and planemos are 
expected. 

In this regime we have 250 objects out of which 103 
have 16 < J < 19 and are thus sufficiently faint to be 
planemos (J = 16 corresponds to ~ 0.015 M Q for zero 
extinction in COND/DUSTY isochrones) and sufficiently 
bright to be accessible for spectroscopy. About one quar- 
ter of them has been observed spectroscopically. Thus, 
the catalogue could contain a few additional planemos - 
up to 6, statistically, which would yield a total of up to 
8. This analysis shows that one particular bias in the 
selection of candidates, the cutoff at II — 12 = 0.3 which 
was necessitated by the large amounts of contaminating 
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objects at bluer colours, might affect our ability to find 
planemos. 

We should not miss any significant number of planemos 
with enhanced mid-infrared colours due to the presence 
of a disk. If the disk fraction for planemos is as high 
as for brown dwarfs, as show n for the cr Ori cluster 
l|Scholz fe Javawardhanal[2008t ). excess emission should 
be expected for about half of the planemos in NGC1333 
(see SONYC-IV for more information on disk fractions). 
These objects should be located to the right of the dashed 
line. In this area we find 39 objects from which about 
half have been verified spectroscopically (including the 
two planemos). This implies that we might miss about 
two more objects in this colour domain. Assuming a disk 
fraction of around > 50%, this again gives a total number 
of up to 8 planemos. 

The JKS catalogue can further be used to probe the 
regime J > 19 which is not accessible for spectroscopy. 
The JKS catalogue contains 292 objects with J > 19, 124 
in the colour space expected for planemos (dash-dotted 
lines in Fig. 1101) . However, this sample is dominated 
by strongly reddened objects with high J — K colour. 
Only about 30 of them have J — K colours consistent 
with Ay < 5 and only 7 are below Ay = 3. Given that 
two thirds of the confirmed brown dwarfs have an opti- 
cal extinction below 5 mag and about half below 3 mag, 
this basically rules out that the candidates too faint for 
spectroscopy down to J ~ 21 contain a sizable popula- 
tion of planemos. According to the COND and DUSTY 
isochrones, this limit corresponds to object masses of 
0.002-0.003 M Q (for A v of 0-5 and age of lMyr). 

5.4. The frequency of planemos 

In the following we will use the survey results in 
NGC1333 to give a quantitative constraint on the fre- 
quency of planemos which can be compared with other 
star forming regions. To do this, we estimate the fraction 
of brown dwarfs which have a mass in the planemo do- 
main, fp = Np/Nbd- This quantity serves as a proxy for 
the shape of the mass function in the substellar regime. 
Here we treat the population of planemos as a subsample 
of the population of brown dwarfs. 

Our survey in NGC1333 yields only 3 objects with 
effective temperatures below 2500 K or spectral types 
later than M9. According to the currently available 1- 
5Myr isochrones these objects have masses below the 
Deuterium burning limit of 0.015 Mq and are thus good 
candidates for planemos. As outlined above, we might 
miss about 6 more planemos due to our cutoff in the 
mid-infrared colour. 

The total number of brown dwarfs can be calculated 
as follows: Adding the objects listed in SONYC-IV and 
this paper, there are 58 with spectral type M5 or later 
or effective temperatures of 3200 K or below. From these 
around 10-20 are likely to be above the substellar bound- 
ary, which leaves 38-48 brown dwarfs. Following the same 
statistical arguments that we have used in Sect. 4.1. of 
SONYC-IV we estimate that we might miss up to 13 
more brown dwarfs, plus 6 planemos (see above), i.e. the 
total is 57-67. 

Thus, for NGC1333 we estimate that f P is 12-14%. 
Just counting the number of confirmed objects, the ra- 
tio becomes 6-8%. This includes all planemos down to 
masses of 0.006-0.008 Mq (the mass limit for J = 19 and 



A v = 5). 

This value can be compared with other star forming 
regions. This comparison comes with three caveats: 1) 
Only a Ori and UpSco have been surveyed with a depth 
comparable to our survey in NGC1333, i.e. all other 
ratios should be considered lower limits. 2) While the 
general survey strategy is similar in all cases (candidate 
selection based on optical/near-infrared photometry plus 
follow-up spectroscopy), the details of the strategy differ 
from region to region, which might introduce unknown 
biases. To be conservative, we focus in the following 
mostly on samples with spectroscopic confirmation. 3) 
Most of the studies listed below do not make an attempt 
to estimate the total substellar population, by taking into 
account possible survey biases. 

We note that all existing surveys use the 
COND/DUSTY isochrones for mass estimates, as 
in our paper. In the case these isochrones turn out 
to be systematically off, the results of the comparison 
should still be valid and merely shifted to different mass 
regimes. 

a Ori: The summary paper bv lCaballero~et"aH (|2007f) 
lists 46 confirmed brown dwarfs, out of which 12 are ex- 
pected to be in the planemo regime, i.e. fp ~ 26%. Not 
all of these objects are spectroscopically confirmed, but 
the contamination in this cluster is generally low due to 
the negligible extinction. The surveys in this region are 
similarly deep as in NGC1333, i.e. the numbers should 
be comparable. 

Upper Scorpius: Based on photometry and proper mo- 
tions from UKIDSS, the total number of brown dwarfs 
in this region is 68 in the UKIDSS covered area, out 
of which 49 are in the northern area where follow-up 
spectroscopy is available ([Dawson et al.l[20TTh . The cur- 
rently ongoing spectroscopic verification shows that es- 
sentially all these objects a re indeed members of UpSco 
(jLodieu et al.ll2008l 1201 laf ) . Converting the photometry 
to masses, this sample should cover the mass range down 
to 0.01 Mq and contains about ha lf a dozen objects be - 
low 0.015 M Q . A deeper survey bv lLodieu et alJ(|2011bf ) 
does not find new members in the planemo domain. 1 
This yields fp ~ 10%. 

p-Ophiuchus: This region is affected by very high and 
variable extinction, which strongly hampers any sur- 
vey. The existing census is not nearly complete be- 
low 0.03 Mq. So far, there are around 40 confirmed 
objects with s pectr al types later or eq u al M6 from 
Wilki ng et all (I2008h: iGeers et all d20ll: iMarsh et al " 



(2010 a[): iMuzic et all (|2012f >: lAlves de Oliveira et al . 
(|2010L 120121 " Out of these, 8 might be in the planemo 
regime based on spectral type and effective temperature, 
which results in fp ~ 20%. 

Chamaeleon-I: The c ensus by iLuhmanl (|2007f ): 
iLuhman fe Muenchl (|2008h contains 35 objects with spec- 
tral type of M6 or later which are to be considered 
brown dwarfs. 5 of these have estimated masses below 
0.015 M©, i.e. f p ~ 14%. We note that our own deep 
survey work in Cha-I has not r evealed any addi tional 
planemos yet down to 0.008 M (jMuzic et al.ll2011h . 

Taking into account the caveats listed above, the val- 

10 Here we neglect the 5 T-dwarf candidates identified based on 
imaging in Methane-sensitive narrow-band filters without spectro- 
scopic or proper motion confirmation. 
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ues for fp in the various regions are consistent with each 
other and in the range between 10 and 20%; only erOri 
might have a slightly larger planemo fraction. Focus- 
ing on the confirmed objects, NGC1333 has a relatively 
low planemo fraction compared with other star forming 
areas. 

We can also estimate the fraction of planemos with re- 
spect to the population of young stars, gp = Np/Ns- 
For NGC1333 we have recently determined that the 
star-vs-brown dwarf ratio, defined as R\ = 7V(0.08 — 
1.0 M o /iV(0.03 - 0.08 M Q ), is 2.3 ± 0.5, for most other 
regions this ratio is larger. Taking this into account, for 
NGC1333 the fraction of planemos with respect to the 
total cluster population is in the range of 2-3%. For 
the other regions where this type of census is available 
this number becomes 2-5%, i.e. stars are 20-50 times 
more frequent than planemos. This is consistent with 
the pla nemo fraction of < 10% ()Lucas et alj|2005l) and 
1-14% (jLucas et alJl2006t ) derived for the Orion Nebulae 
Cluster (ONC), albeit for slightly different mass limits. 
Thus, the number of planemos is insignificant compared 
with the total population of stars and brown dwarfs. The 
planemo contribution to the mass budget of young clus- 
ters is negligible (well below 1%). 

5.5. The mass spectrum 

With the available database of spectroscopically con- 
firmed objects in NGC1333, we are in the position to 
construct the Initial Mass Function (IMF) in the low- 
mass regime. We choose to express the IMF as a mass 
spectrum dN/dM oc M~ a . To do that, we use the same 
spectroscopic samples as in S ect. 15.11 combining th e con- 
firmed low-mass stars from IWinston et al.l (|2009T) with 
the census of very low mass objects. We select all M and 
L dwarfs, in total 107 objects, out of which 100 have an 
estimate for t he effective temperatu re, either from our 
survey or from IWinston et al.l (|2009H . 

Masses were derived for this sample by comparing the 
temperatures with theoretical evolutionary tracks. Tem- 
peratures do not significantly depend on age in the 1- 
10 Myr phase, as opposed to luminosities and magni- 
tudes, and thus provide a mass estimate that is insensi- 
tive to age spread. We use the 1 Myr isochrone from the 
BCAH models for M > 0.1 Mq and the 1 Myr isochrone 
from the DUSTY models for M < 0.1 M©. For each 
object with a given T e g, we selected the isochrone data- 
points within ±200 K around T e g, fitted this section ei- 
ther linearly (if there are less than 5 datapoints) or with 
a second order polynomial, and used that function to 
calculate the mass corresponding to T e g. All masses are 
tied to this isochrone and should not be compared with 
values derived using other models. 

The mass spectrum was calculated in a way to achieve 
similar statistical uncertainties in each mass bin, i.e. with 
varying bin size. For M > 0.7 Mq the sample sizes are 
too small for a meaningful analysis, the 10 objects above 
this limit were excluded. The lowest mass bin includes 
all objects with M < 0.03 Mq. In Fig. [TJJwe plot our 
result. The datapoints in this diagram are consistent 
with dN/dM oc M~ a with a = 0.61, overplottcd as 
dashed line. For the substellar regime alone, the best- 
fitting slope is a = 0.46. For the stellar regime above 
0.1 Mq the slope becomes a = 1.0. The uncertainty in 
the slope is in the range of 0.1. 
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Fig. 11. — Mass function in NGC1333 based on the spectroscop- 
ically confirmed members with measured effective temperatures. 
The binsizes were chosen to achieve approximately equal error 
bars. The solid lines show the Kroupa mass function (a = 1.3 
for 0.08 < M < 0.5M Q , a = 0.3 for M < O.O8M ). The dashed 
line shows the best fit to the data (a = 0.6). 

To verify how the uncertainties in the effective temper- 
atures (±200 K) affect these results, we carried out the 
following test: For half the objects, randomly selected, 
we increased T c g by 200 K, for the other half we reduced 
it by 200 K. This effectively scrambles the values within 
the errorbars. Then the mass spectrum was calculated 
again using the same procedure. This yields a = 0.66 for 
the entire sample and 0.55 for the brown dwarfs, slightly 
higher as before but well within the errorbars. This test 
also shows that the peak around 0.2 Mq seen in the figure 
is spurious. 

Also shown is the Kroupa IMF, an empirical descrip- 
tion of the mas s spectrum as a segmented power law, 
as a solid line (jKroupal l2001| ) . This function follows 
a = 1.3 ± 0.5 in the stellar regime up to 0.5 M and 
a = 0.3 ± 0.7 in the substellar domain. Within the un- 
certainties, the Kroupa IMF is consistent with the mass 
spectrum in NGC1333. 

The slope derived here is broadly consistent with val- 
ues derived for most other star forming regions and young 
clusters. Some examples: For the very low mass objects 
in a Ori, values of 0.6-0.7 have been published i n the lit- 
erature (jCaballero et al. 2007: B eiar et al.ll2011|) . For the 
domain from 0.03 to 0.6 Mq the cl uster Blanco 1 has a 
mass function with a = 0.69 ±0.15 (jMoraux et al.ll2007f ). 
For Upper Scorpius and the mass regime .01-0.3 Mq, a 
value ofa = 0.6±0.1 has been published (|Lodieu et al.l 
I2007h . For the Pleiades a is found to b e in the range 
of 0. 6-1.0 in the low-mass regime (e.g. iMoraux et al.1 
120031 and referenc es therein) . For the a Per cluster , the 
slope is 0.59±0.05 (|Barrado v Navascues et aTll2002D . In 
Collinder69 the slope has been f ound to be betwe en 0.2 
and 0.4 for the low- mass regime (|Bavo et al.ll201lf) . 

To our knowledge, the only young region where a sig- 
nificantly different slope of the mass spectrum has been 
determ ined is the O r ion N ebulae Cluster. In a recent 
paper, iDa Rio et ail (|201lD find T of -1.1 to -3.1 for 
the substellar regime, which corresponds to a negative 
a (a = r + 1). Similar results have been found i n 
previous studies in the ONC (e.g. iMuench et al.l [2002). 
These numbers appear to be anomalous compared with 
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the other regions mentioned above. 

Apart from the overall slope, we can examine the shape 
of the mass spectrum in NGC1333. In equal sized bins, 
the numbers per bin increase continuously with decreas- 
ing mass, with the exception of a slight minimum around 
0.1-0.2 Mq. However, the test with the modified temper- 
atures (see above) indicates that this minimum is spuri- 
ous, the data is indeed consistent with a monotonically 
rising histogram. The same result has been obtained 
for Upper Scorpius and a Ori, while some other regions 
show a peak in the mass spectrum around 0.1-0.2 M Q 
and a smaller slope for lower masses (IC348, Ch a-I, p- 
Oph. lLuhman1l2007t lArves de Oliveira et alJl2012f ). One 
consequence of this difference is a smaller star-to-brown- 
dwarf ratio (i.e. a large number of brown dwarfs relative 
to stars) in regions with monotonically rising mass spec- 
trum - as it is indeed found for NGC1333. As discussed 
in SONYC-IV, a comparison of these ratios possibly indi- 
cates environmental differences in the formation of brown 
dwarfs. 

The focus of this study is the planetary mass domain. 
Assuming a monotonic continuation of the power law 
shown in Fig. QTJwith a = 0.6±0.1, we would expect 8±3 
objects with 0.005 < M < 0.015 M , i.e. in the planemo 
domain and detectable in our survey. For comparison, 
in Sect. 15.31 we derive an upper limit of 8 planemos in 
this mass regime from our survey, taking into account the 
incompleteness in the spectroscopic follow-up of our can- 
didates. Thus, the data is consistent with a monotonic 
power law mass spectrum across the Deuterium burn- 
ing limit, but a smaller slope is possible as well. Our 
survey rules out that the slope of the mass spectrum in 
NGC1333 increases in the planemo regime. 

Again this is in agreement with previous findings 
for other regions. In a Ori the census of sub- 
stellar objects is well approxim ated by a monotonic 
slope down to p l anetary masses (ICaballero et all [20071: 
iBeiar et "aLll2011fc iPefia Ramirez et al.ll2012D . with a pos- 
sible turnover (i.e. a lower a) in the mass spectrum 
below 0.006 Mq. F or Upper Scorpius, the results by 
lLodieu et al.l (120 lib!) favour a tur ndown below 0.01 M Q . 
We note that iMarsh et al.l (|2010bt) show a mass function 
of p-Ophiuchus covering the entire planemo regime based 
on a photometrically selected candidate sample, which 
resembles the one in a Ori. Taken together, it seems safe 
to conclude that the slope of the mass function in the 
planemo domain is a < 0.6. 

We note that a < 0.6 corresponds to a planemo frac- 
tion as defined in Sect. 15.41 of fp < 20% (if the entire 
planemo domain is considered) or fp < 30% (with an 
lower mass cutoff at 0.005 Mq). Indeed, no regions shows 
a larger planemo fraction, most surveys indicate signifi- 
cantly lower values, including our own work in NGC1333. 

Thus, the current census of star forming regions pro- 
vides support for the idea that planemos are an exten- 
sion of the population of stars and brown dwarfs and 
form through the same mechanisms. The observations 
indicate that free-floating objects with planetary masses 
exist down to at least 0.006 M Q , but, as pointed out in 
Sect. 15 .41 their numbers and their mass budget are in- 
significant compared with the total population of stars 
and brown dwarfs. 

Some theoretical estimates for the lower mass limit 
in the star formation process - the opacity limit for 



fragmentation - are .007 (|Low fc Lvnden-Belll 11976ft 
or 0.01 (ISilld I1977T) from analytical argument s and 

< 0.001 (IBossII2001D o7 0.003-0.009 M Q (lBatdl2005ft from 
numerical simulations. Except for the lowest-mass val- 
ues, all these estimates are consistent with the current 
o bservational pictur e. 

ISumi et ah! (|2011[ ) recently determined the slope of 
the mass spectrum based on microlensing surveys, and 
find 1.3±g;| for M < 0.01 M©, significantly higher than 
for brown dwarfs (0.49 in their work), indicating that 
Jupiter-mass objects are almost twice as common as 
stars. Their estimate includes planetary-mass objects 
without host stars (i.e. free-floating) and those on wide 
orbits, but they argue that the majority o f them are in- 
deed unbound (but see lQuanz et al.l ()2012l ) for a criticism 
of this claim) . This is in clear disagreement with the re- 
sults from the surveys in star forming regions, where the 
number of planemos per star is much smaller (see Sect. 
I5.4|) and the slope in the mass spectrum is not rising 
(see above). If the microlensing result holds, the addi- 
tional planemos would have to be in the mass regime 

< 0.005 Mq which is not sufficiently probed by the clus- 
ter surveys. This would imply a break in the slope of 
the mass spectrum in the planemo domain, favouring a 
scenario where the Jupiter-mass objects form in a dif- 
ferent way than the planemos and brown dwarfs seen 
in the star formi ng regions - as indeed concluded by 
ISumi et ail (|2011l) . Larger samples of short-period mi- 
crolensing events and deeper surveys in star forming re- 
gions are critical to understand if they are in fact different 
populations or not. 

In another recent paper. IStrigari et ah! (|2012l ) estimate 
that there may be up to 10 5 objects with masses be- 
tween 10 -8 and 10 -2 M Q (called 'nomads' in their pa- 
per) per main-sequence star. To arrive at this number, 
they use the slope of the mass spectrum in the brown 
dwarf/ planemo regime as a constraint. Their upper limit 
of 10 can only be obtained by using a — 2. As shown 
above, this is inconsistent with the surveys in star form- 
ing regions. For a more realistic value of a = 0.5, the 
result given in their paper is > 1 nomads per main- 
sequence star. Given the enormous difficulties in extrap- 
olating the mass function over 6 orders of magnitude in 
mass and different modes of formation, these estimates 
remain highly speculative. 

6. CONCLUSIONS 

We present new spectroscopy for 19 very faint can- 
didate members in the young cluster NGC1333. 7 of 
them are confirmed as very low mass objects, most or 
all of them should be brown dwarfs. Based on the com- 
bined spectroscopic follow up from SONYC-I, SONYC- 
IV and this paper, we are able to put limits on the IMF 
and the number of free-floating planetary-mass objects 
(planemos, M < 0.015 Mq) in this cluster. The mass 
spectrum is in good agreement with a monotonically ris- 
ing power law dN/dM oc M~ a , with a — 0.6 ± 0.1 for 
M < 0.6M e . The overall slope is in agreement with all 
previous studies in other regions, with the ONC as only 
exception. 

Among the substellar population in NGC1333 are 3 
with estimated masses in the planemo domain. Tak- 
ing into account the incompleteness of the spectroscopic 
follow-up in our survey, we estimate that the cluster 
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could harbour a total of 8 planemos down to masses of 
~ 0.005 Mq. This translates into a planemo fraction (the 
fraction of brown dwarfs in the planetary mass domain) 
of 12-14%. This is comparable to the planemo fractions 
determined for other regions. We estimate that there are 
20-50 times more stars than planemos in young clusters. 
Based on the constraints on the number of planemos in 
star forming regions, the slope of the mass spectrum is 
< 0.6 below 0.015 M , i.e. it does not increase in the 
planemo domain. These findings support the idea that 
planemos form through the same mechanisms as young 
stars and brown dwarfs. 
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